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基于 噪声 相关 性 的 惩罚 加 权 最 小 二 乘 算 法 在 低 剂 量 数字 乳腺 层 


析 成 像 中 的 应 用 
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摘要 :目的 对 投影 数据 方差 的 精确 建 模 并 结合 DBT 平 板 探测 系统 的 噪声 相关 性 构建 精准 噪声 模型 下 的 基于 噪声 相关 性 的 惩 
罚 加 权 最 小 二 乘 算法 在 低 剂量 乳腺 层 析 成 像 图 像 中 的 应 用 。 方 法 首先 对 投影 数据 进行 量子 噪声 和 电子 噪声 建 模 , 使 以 往常 用 
的 近似 噪声 模型 精准 化 ,然后 构建 基于 噪声 相关 性 的 惩罚 加 权 最 小 二 乘 算法 用 于 投影 数据 恢复 ;最 后 对 处 理 后 的 投影 数据 采用 
滤波 反 投 影 算 法 进行 重建 。 结 果 对 不 同 剂 量 下 ACR 标 准 体 模 数 据 进行 处 理 得 到 的 重建 结果 噪声 明显 降低 ,细节 对 比 度 提 


高 。 恢 复 投影 
显 ,重建 DBT 图 像 质量 有 很 大 的 提升 。 


关键 词 :数字 乳腺 层 析 成 像 ; 低 剂量 ;噪声 相关 性 ;加 权 最 小 二 乘 


据 的 重建 图 像 与 原始 数据 重建 图 像 相 比 ,CNRs 和 LSNRs 提 升 了 约 3.6 倍 。 结 论 对 投影 数据 噪声 抑制 效果 明 
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Abstract: Objective To achieve low-dose digital breast tomosynthsis (DBT) projection recovery using penalized weighted least 


square algorithm incorporating accurate modeling of the variance of the projection data and noise correlation in the flat panel 
detector. Methods Models were established for the quantal noise and electronic noise in the DBT system to construct the 


penalized weighted least squares algorithm based on noise correlation for projection data restoration. The filter back projection 
algorithm was then used for DBT image reconstruction. Results The reconstruction results of the ACR phantom data at 
different dose levels showed a good performance of the proposed method in noise suppression and detail preservation. CNRs 


and LSNRs of the reconstructed images from the restored projections were increased by about 3.6 times compared to those of 
reconstructed images from the original projections. Conclusion The proposed method can significantly reduce noise and 


improve the quality of DBT images. 
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近年 来 ,乳腺 瘤 的 发 病 率 一 直 呈 现 上 升 趋势 ” , 乳 
腺 筛 查 成 为 了 降低 乳腺 癌 死 亡 率 的 重要 手段 ” 。 

义 线 乳腺 摄影 技术 “对 微小 病灶 具有 高 度 敏感 性 ， 
但 传统 的 全 视野 数字 乳腺 XX 摄影 (FFDM) 由 于 平面 二 
维 成 像 的 组 织 重 受 效应 在 乳腺 肿瘤 的 检测 中 存在 一 定 
程度 的 检 出 失败 率 *9, 同 时 ,由 正常 组 织 重 闭 而 产生 的 
假 阳性 肿瘤 会 误导 医生 的 诊断 。 近 年 来 ,数字 乳腺 层 析 
成 像 (DBT) 逐 渐 在 临床 中 得 到 广泛 应 用 区 ,DBT 成 像 


利用 传统 的 X 线 源 , 以 平板 探测 器 为 旋转 中 心 采 集 一 定 
数量 的 投影 数据 ,并 通过 有 限 数 量 的 投影 数据 进行 三 维 
重建 ,实现 了 乳腺 的 三 维 成 像 。 数 字 乳 腺 层 析 成 像 很 大 
程度 上 解决 了 二 维 成 像 组 织 重 二 所 导致 的 检 出 失败 和 
假 阳 性 误诊 的 问题 ”。 同 时 ,相对 于 CT 成 像 技 术 ,DBT 
成 像 很 大 程度 上 降低 了 患者 的 辐射 剂量 ,但 乳腺 作为 
一 个 对 X 射 线 敏 感 融 官 ", 其 辐射 危害 仍然 不 容 小 凯 。 
为 进一步 降低 辐射 剂量 ,可 以 通过 降低 X 射 线 球 管 电流 
(mAs) 来 实现 超 低 剂 量 DBT 成 像 中 ,但 超 低 剂量 扫描 条 
件 下 的 重建 图 像 噪声 水 平 较 高 ,导致 低 对 比 度 病灶 及 微 
小 钙化 点 难以 识别 。 

为 抑制 DBT 图 像 噪声 ,提高 病灶 检 出 率 ,本 文 借鉴 
低 剂量 CT 成 像 中 的 投影 数据 降 噪 理论 “对 投影 数据 
进行 恢复 后 重建 ,并 在 噪声 模型 构建 和 恢复 算法 设计 两 
个 方面 做 了 改进 ,以 此 提升 投影 数据 恢复 效果 和 重建 图 
像 质量 。 在 噪声 模型 构建 方面 ,Ma 等 "在 对 投影 数据 
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的 方差 分 析 过 程 中 发 现 低 剂量 扫 描 条 件 下 电子 噪声 的 
忽略 会 使 测量 投影 数据 的 方差 估计 与 真实 值 产生 较 大 
差距 ,并 提出 和 验证 了 精确 化 的 噪声 模型 ,本 文 基于 这 
一 理论 将 常规 的 近似 噪声 模型 做 了 精确 化 处 理 , 用 于 低 
剂量 DBT 投 影 数据 恢复 。 同 时 ,我 们 发 现 基于 平板 探 
测 需 获取 的 投影 数据 存在 噪声 相关 性 “"。 基 于 上 述 ,本 
文 充 分 利用 DBT 平 板 探测 器 的 数据 统计 特性 ,在 对 
DBT 投 影 数据 进行 精确 噪声 建 模 ( 对 投影 数据 进行 方 
差 特性 统计 ) 的 基础 上 ,结合 噪声 相关 性 构建 惩罚 加 权 
最 小 二 乘 算法 对 投影 图 像 进行 数据 优化 。 最 后 利用 改 
进 的 FBP 算 法 进行 图 像 重 建 ,以 提升 DBT 图 像 质 量 。 


1 方法 
1.1 投影 数据 噪声 建 模 

对 于 图 像 噪 声 的 准确 建 模 ,探测 器 响应 的 数据 统计 
尤为 重要 。Whiting 等 人 对 能 量 积分 型 探测 器 数据 校正 
以 及 数据 采集 的 随机 过 程 进行 理论 建 模 ,得 出 校正 后 数 
据 满足 复合 泊 松 分 布 的 结论 "2 ,透射 光子 分 布 函数 将 
量子 噪声 与 背景 电子 噪声 结合 ,定义 为 : 

T, ~ X G(E,)Pois son(1(E,)) + Normal{d o.) (1) 


i 表示 探测 单元 序列 编号 BR k RIR KICAEEK , 
GO 表示 能 谱 概 率 , Possion{.} 表征 数据 的 泊 松 统计 特 
性 ,其 中 4 代表 穿 过 测量 体 后 的 光子 平均 数 , 即 探测 器 
实际 接收 到 的 信号 。 公 式 (1) 的 第 二 项 表征 电子 品 声 ， 
服从 均值 为 d, ,方差 为 o.r 的 高 斯 分 布 ,电子 噪声 的 参 
数 可 以 通过 暗 电流 数据 测量 。 对 于 复合 泊 松 分 布 ,目前 
还 还 缺 闭合 准确 的 似 然 函 数 ,本 文 忽略 久光 多 能 特性 ， 
将 X 线 光子 分 布 函数 简化 为 : 

I, ~ Pois son(4) + Normal (d, c, P} (2) 

利用 朗 伯 比尔 定律 对 探测 需 数 据 进行 对 数 变换 ,可 
以 得 到 探测 单元 i 对 应 路 径 上 的 线性 衰减 系数 线 积 人 
数据 : 

-学 -km7 (3) 


ROP, L WAWA £s UE I B ETA 
目 ,可 以 通过 空 扫 描 进行 测量 ; 7 表示 探测 器 接收 到 的 
透射 光子 量 的 随机 变量 ,分 布 函 数 如 公式 (2) 所 示 。 忽 
略 元 的 随机 特性 , p, 是 一 个 随机 变量 ,其 均值 方差 可 
以 通过 计算 求 取 。 本 文 线 积分 投影 数据 方差 模型 ”: 

E 1 0 -1.25 

a =var(p)= " (L+ i ) (4) 

IF 4,21, expo p) 。 式 (4) 给 出 了 投影 数据 方差 
和 光子 平均 数 .电子 噪声 方差 的 关系 ,保留 了 量子 噪声 
和 电子 噪声 对 方差 特性 的 影响 。 
1.2 基于 嗓 声 相关 性 的 惩罚 加 权 最 小 二 乘 算法 在 DBT 
投影 数据 恢复 中 的 应 用 

基于 先前 对 平板 探测 需 的 噪声 相关 性 的 研究 ,在 同 


一 扫描 条 件 下 ,对 大 量 重复 的 平板 探测 顺 投 影 数 据 统计 
分 析 , 得 出 相 邻 投影 数据 之 间 存 在 着 噪声 相关 性 的 结论 ， 
不 同 路 径 积分 数 的 噪声 相关 性 可 由 相关 因子 B, 描述 : 

av。 cov(p,p) m 

|99, vap) var(p) 

基于 式 (5) 中 测定 的 相关 系数 ,利用 式 (4) 中 对 投影 
数据 噪声 方差 的 估计 ,我 们 可 以 获得 
cov(p.p) c, - p, |var(p) |var(p) ,由 此 构建 投影 数 
据 噪 声 协 方差 矩阵 Y, ,该 矩阵 是 带宽 为 8B+1 实 对 称 对 
角 和 矩阵 , 即 : 当 isj, i-j- 2l, i-j-£B, i-j- «B4 1) 
时 ,矩阵 > 的 元 素 o; 为 非 零 值 ,矩阵 的 其 他 元 素 均 取 值 
为 零 , 当 采取 逐 列 方式 重 排 数 据 时 ,B 等 于 列 维 度 , 当 采 
用 逐 行 重 排 时 ,8 等 于 行 维度 。 基 于 噪声 相关 特性 , 利 
用 加 权 最 小 二 乘 算法 (PWLS) 以 投影 数据 建立 能 量 函 
数 ,估计 最 佳 线 积 4 

dXp) -(? - p? G - p) * BR) (6) 

式 46) 中 第 一 项 为 数据 保 真 项 ,其 中 了 表示 探测 天 
直接 测量 值 对 数 变换 后 的 数据 , p 表示 待 估计 的 最 佳 线 
积分 值 , > 为 线 积 分 投影 数据 的 协 方差 矩阵 。 公 式 (6) 
中 第 二 项 为 惩罚 项 ,其 中 p 为 调节 保 真 项 和 签 罚 项 比例 
的 参数 ,惩罚 项 设置 如 下 : 

R(p)= Rp)= YY wp, -pe) (7) 


S, 指 在 二 维 投影 数据 中 与 像素 /最 相 邻 的 4 个 像 
素 , 本 文中 4 个 像素 的 权 w(hj) 取 相 同 值 。 通 过 最 小 化 
代价 函数 (6) 可 以 得 到 优化 的 线 积分 投影 数据 : 

p argmin($-pJE G-p) + 

BY. we, p) 


本 文中 协 方差 矩阵 的 非 负 对 称 特性 以 及 惩罚 先 验 
项 R(p) 的 二 次 平滑 设计 使 得 式 (8) 中 的 代价 函数 为 辆 
函数 。 由 凸 函数 的 性 质 , 我 们 可 以 将 对 B(p) 的 最 优化 
转化 为 求解 VB(p)=0 的 极 值 问题 ,上 且 该 解 理论 上 是 唯 
一 的 。 将 VR(p) P828 VR(Q) - V|p|| ^ Ap. A29" 
JE Zr EAE KESZ DAS TOSE FA RERE s AE : 


(8) 


4, i=j 
A;=53-l, i-j=lori-j=+B (9) 
0, otherwise 


因此 代价 函数 的 梯度 形式 可 以 通过 下 式 改 写 : 

Vo p)- V -p) X G-p) * ENRQ)-0 

X (p-y)*fAp-0 

FRU, RATI LORIE 7 (p -y)+B4p=0 最 优化 求 
解 式 (8)。 

为 了 充分 利用 和 矩阵 > A A BA iiL ARRA ,我 们 可 以 
利用 LSQR 算 法 进行 迭代 求解 。 对 Y (p-y) -BAp-0 
迭代 求解 过 程 如 下 : 

0. 初 始 化 p^ 


(10) 


JEBOEXSIIR IET € ; 
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1. 求 解 线性 方程 b=y-p”, 以 更 新 b; 
2. 对 步 又 1 求解 得 到 的 b 值 ,求解 线性 方 


TE p4p" =b ,得 到 p"; 

3. 当 满足 p" ^ -p^| «e 停止 迁 代 ,否则 返 
回 步骤 1; 

* 其 中 为 中 间 变 量 


3 恢复 投影 数据 重建 
对 恢复 后 的 投影 数据 ,本 文 利用 三 维 的 滤波 反 投 影 
(FBP) 算 法 进行 重建 ””“。 重 建 算法 基于 滤波 反 投影 ,为 
实现 乳腺 的 最 佳 成 像 效果 改进 滤波 设计 ,同时 采用 距离 
坐标 变换 后 再 进行 反 投 影 (距离 驱动 反 投 影 ,DDBP)”， 


这 样 的 方式 相 比 传统 的 对 整个 图 像 基于 路 径 的 反 投 影 
缩短 了 计算 时 间 , 从 精度 的 角度 出 发 ,距离 坐标 变换 反 
投影 考虑 了 图 像 像素 和 探测 噩 的 宽度 问题 ,可 以 忽略 将 
投影 数据 的 锥 形 束 近似 平行 束 处 理 时 对 重建 的 影响 。 


2 结果 
2.1 实验 数据 获取 

本 文 主 要 采用 DBT 实 验 原 型 机 采集 的 ACR 标准 
体 模 投影 数据 进行 图 像 重 建 。ACR 标准 体 模 是 乳腺 成 
像 所 使 用 的 标准 测试 体 模 , 体 模 内 部 包含 与 乳腺 结构 相 
似 衰 减 系数 的 物质 ,对 应 于 乳腺 的 腺 体 .纤维 组 织 .钙化 
点 `. 圳 肿 等 结构 ,ACR 标准 体 模 结 构 信 息 参 见 图 1。 


图 1 乳腺 成 像 标准 体 膜 ACR 体 模 内 部 结构 信息 及 实物 图 


Fig.1 Internal structure and appearance of an ACR standard phantom. 


数据 采集 系统 为 Analogic 公 司 AXS-2430 平 板 探 
测 器 ,投影 数据 的 采集 参数 为 : 管 电压 30 kV , 管 电流 分 
别 是 40.60.80.100 mA。 探 测 絮 以 与 X 线 源 曝光 相同 
的 时 序 采集 25 个 投影 ,X 线 源 在 +24° 范 围 内 每 隔 2° 曝 光 
1 次。 数据 采集 的 过 程 中 ,XX 线 源 到 探测 器 的 距离 为 
650 mm, 旋 转 中 心 到 探测 器 的 距离 为 25 mm。 原 始 投 
影 数 据 大 小 为 3584x2816, 像 素 大 小 为 0.085 mmx 
0.085 mm, 

22 评价 参数 

为 定量 评价 算法 性 能 ,本 文中 我 们 选用 局 部 信 噪 比 
和 对 比 度 信 品 比 作为 图 像 质量 的 评估 指标 ,其 中 ,局 部 
信 噪 比 LLSNR) 和 对 比 度 信 噪 比 (CNRJ) 分 别 定义 为 : 


LSNR - flo, (11) 
2| f-b 

CNR- TL (12) 
0, t 0, 


和 背景 区 域内 的 像素 均值 , o, , o, 为 感 兴趣 区 域 (ROD 
和 背景 区 域内 的 标准 差 。 
2.3 投影 图 像 恢复 前 后 比较 

本 文 的 方法 实现 过 程 中 ,首先 通过 暗 扫描 对 电子 噪 


声 方 差 进行 估计 ,再 结合 公式 (4) 进 行 投影 数据 的 方差 
估计。 图 2 给 出 了 0 角度 下 利用 基于 噪声 相关 性 的 惩罚 
加 权 最 小 二 乘法 恢复 前 后 的 投影 图 像 。 其 中 图 2 PAB 
的 扫描 参数 为 30 kV .40 mA, C.D 的 扫描 参数 为 30 kV、 
60 mA, KIA 和 图 C 为 被 噪声 污染 的 投影 图 像 ,图 B 和 
图 D 为 利用 本 文 算法 恢复 后 网 像 。 

对 比 不 同 剂量 下 的 投影 图 像 ,可 以 看 出 低 剂 量 扫描 
会 导致 投影 数据 存在 严重 的 噪声 ,利用 基于 噪声 相关 性 
的 惩罚 加 权 最 小 二 乘 算法 , 低 剂量 投影 图 像 的 噪声 得 到 
有 效 的 抑制 ,投影 图 像 无 失真 。 

图 3 展示 了 30 kV,40 mA 剂量 下 的 0 角度 投影 
据 在 参数 p= 1200 时 本 文 算 法 迭代 的 收敛 情况 , 横 轴 表 
示 迁 代步 数 , 纵 轴 量 | -由 可 以 表征 第 二 步 与 第 
n-1 步 迭代 后 投影 数据 结果 之 间 的 差异 ,从 图 中 可 以 看 出 
本 文 算法 收敛 速度 较 快 ,10 步 迭代 即 可 获取 稳定 结果 。 

本 文中 ,我 们 对 0 角度 下 的 投影 数据 进行 了 LSNR 
对 比分 析 , 图 4 中 折线 图 展示 了 不 同 剂量 投影 图 像 恢复 
前 后 的 LSNR 值 ,从 图 中 我 们 可 以 看 出 ,基于 噪声 相关 
性 的 惩罚 加 权 最 小 二 乘 算 法 可 以 有 效 地 抑制 DBT 投 影 
图 像 噪声 ,提高 投影 图 像 的 信 品 比 。 
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2 30 kVp、40 mA 和 30 kVp、60 mA 剂量 下 0 角度 投影 数据 恢复 前 后 的 图 像 
Fig.2 Original and restored projection images at angle 0 with different scan settings. A: Original projection image obtained at 30 
kVp and 40 mA; B: Restored projection image obtained at 30 kV and 40 mA; C: Original projection image obtained at 30 kV and 


60 mA; D: Restored projection image obtained at 30 kV and 60 mA. 
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Fig.3 Convergence of PWLS algorithm at 30 kVp and 40 
mA at the angle of 0 with B=1200. 


2.4 重建 结果 比较 

利用 距离 驱动 反 投影 (DDBP)" 算 法 对 恢复 后 投影 
图 像 进行 重建 ,我 们 可 以 得 到 如 图 5 所 示 的 重建 图 像 。 
图 5A~ 五 分 别 给 出 了 40、60、80、100 mA 四 个 剂量 下 原 
始 投 影 数 据 重建 图 像 和 恢复 投影 数据 重建 图 像 。 

从 重建 结果 中 我 们 可 以 看 出 ,利用 原始 投影 数据 重 


别 ,同时 图 像 的 对 比 度 也 十 分 有 限 ;而 利用 恢复 投影 数 
据 重 建 的 图 像 中 ,噪声 水 平 较 低 ,点 状 细节 得 到 凸显 。 
相对 于 图 A 和 C, 图 B 和 D 中 第 三 级 点 状 结构 可 以 很 好 
的 重建 并 显现 出 来 ,相对 于 图 E 和 G, 网 F 和 五 中 第 四 
级 点 状 结构 可 以 很 好 的 重建 并 显现 出 来 ;由 此 可 见 ， 
利用 基于 噪声 相关 性 的 惩罚 加 权 最 小 二 乘 算法 可 以 实 
现 低 剂量 DBT 图像 的 有 效 去 噪 , 提 高 DBT 图 像 的 对 
比 度 。 

为 展示 本 文 算法 的 有 效 性 ,我们 采用 BM3D 
(Block Matching 3-Dimension ) 滤 波 算法 只 2 用 于 投影 
数据 恢复 ,然后 对 恢复 数据 进行 FBP 重 建 。 图 6 中 展示 
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图 4 0 角度 下 40、60、80、100 mA 四 个 剂量 的 投影 数据 恢 
复 前 后 局 部 信 噪 比 折线 图 

Fig.4 Local SNRs calculated from projection data at angle 
0 and at the dose levels of 40, 60, 80 and 100 mA before 
and after denoising. 


T 40 mA 剂量 FFBM3D 及 本 文 方法 重建 的 DBT 图 像 ， 
可 以 看 出 图 6B 中 点 状 结构 清晰 度 明 显 优 于 图 6A 对 应 
级 数 点 状 结构 。 此 实验 结果 说 明 , 本 文 的 算法 在 图 像 细 
节 结 构 保 持 方面 有 上 佳 表 现 。 

为 进一步 说 明 方法 的 有 效 性 ,我 们 分 别 计算 了 四 个 
剂量 下 本 文 算法 重建 图 像 与 图 7 所 示 的 相同 感 兴趣 区 
的 LSNRs 和 CNRs。 表 1 和 表 2 分 别 给 出 了 不 同 剂量 下 
投影 数据 恢复 前 后 的 重建 图 像 的 LSNR 及 CNR 值 。 从 
表格 中 数据 可 以 看 出 , 随 着 剂量 的 升 高 ,重建 DBT 图 像 
的 LSNR 及 CNR 值 变 高 ;相对 于 原始 投影 数据 重建 结 
果 , 利 用 PWLS 算 法 恢复 投影 数据 重建 图 像 有 更 高 的 
LSNR 和 CNR。 


3 讨论 

乳腺 筛 查 的 重大 需求 决定 了 数字 乳腺 层 析 成 像 这 
一 新 兴 义 射线 成 像 技 术 在 临床 应 用 上 的 重要 意义 。 然 
而 ,噪声 的 存在 可 能 掩盖 某 些 低 对 比 度 的 病灶 区 域 ,在 
图 像 断 层 重建 的 去 噪 人 研究 中 ,基于 最 小 二 乘 算法 的 去 噪 
方法 得 到 了 广泛 关注 “”。 本 方法 借鉴 了 前 人 成 果 低 
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5 不 同 剂量 下 投影 数据 恢复 前 后 的 重建 结果 


Fig.5 Reconstruction results before and after projection data recovery at different doses. A: The image reconstructed from 
original projections at 30 kVp and 40 mA; B: The image reconstructed from the restored projections at 30 kVp and 40 mA; C: 
The image reconstructed from original projections at 30 kVp and 60 mA; D: The image reconstructed from the restored 
projections at 30 kVp and 60 mA; E: The image reconstructed from original projections at 30 kVp and 80 mA; F: The image 
reconstructed from the restored projections at 30 KVp and 80 mA; G: The image reconstructed from original projections at 30 
kVp and 100 mA; H: The image reconstructed from the restored projections at 30 kVp and 100 mA. 


图 6 不 同 算法 对 40 mA 投影 数据 处 理 后 的 重建 图 像 . 
Fig.6 Reconstructed images from restored Projection data (40 mA) using different 
algorithms. A: BM3D; B: the proposed method. 


剂量 CT 成 像 中 的 投影 数据 降 噪 算法 ,但 采用 了 更 加 精 
确 地 噪声 模型 并 结合 噪声 相关 性 构建 惩罚 加 权 最 小 二 
乘 算法 对 投影 图 像 进行 数据 优化 ,对 于 噪声 模型 的 改 
进 , Liang 等 人 在 对 投影 数据 的 方差 分 析 过 程 中 发 现 ， 
噪声 分 布 模型 的 正 态 拟 合 与 实际 情况 存在 一 定 差距 ， 
在 扫描 剂量 较 高 的 情况 下 ,探测 器 系统 所 导致 的 电子 哄 
声 对 探测 絮 测 量 的 投影 数据 影响 极 低 , 因 此 可 以 忽略 不 
计 , 但 是 在 低 剂 量 扫 描 条 件 下 ,忽略 电子 噪声 会 使 测量 值 
的 方差 估计 与 真实 值 产 生 较 大 差距 ,本 方法 针对 低 剂量 
DBT 投 影 数据 的 噪声 进行 恢复 ,噪声 模型 的 精确 化 无 
疑 使 降 品 效果 更 加 显著 。 

将 噪声 相关 性 推论 得 到 的 数据 方差 特性 用 于 最 小 
二 乘 算法 中 可 以 优化 数据 恢复 结果 。DBT 成 像 所 用 的 
球 管 及 管 电压 相对 CT 成 像 略 有 不 同 ” ,但 单一 角度 下 


DBT 投 影 数 据 与 CT 投影 数据 的 产生 机 制 并 无 本 质 差 
异 。 在 数据 采集 过 程 中 ,DBT 采 用 了 平板 探测 器 ,不 同 
探测 器 单元 上 的 投影 数据 可 能 相互 影响 。 本 文 对 投影 
数据 进行 方差 特性 统计 ,并 结合 噪声 相关 性 构建 惩罚 加 
权 最 小 二 乘 算法 对 投影 图 像 进行 数据 优化 ,利用 改进 的 
FBP 算 法 重建 ,以 提升 DBT 重 建 图 像 质量 。 人 参照 CT 成 
像 中 噪声 模型 构建 方法 ,将 电子 噪声 和 量子 噪声 结合 构 
建 加 权 最 小 二 乘 算法 来 对 噪声 进行 去 除 。 通 过 合理 的 
代价 函数 优化 及 选取 合适 的 参数 ,投影 图 像 恢 复 结果 及 
重建 DBT 结 果 均 显示 ,基于 噪声 相关 的 惩罚 加 权 最 小 
二 乘 算法 在 低 剂 量 DBT 成 像 中 有 着 上 佳 表现 ,不 同 剂 
量 图 像 的 局 部 信 噪 比 和 对 比 度 信 品 比 提升 约 3.6 倍 。 同 
时 ,为 进一步 说 明 算法 的 有 效 性 ,实验 中 我 们 加 入 了 
BM3D 滤 波 算法 “3 对 投影 数据 进行 去 噪 人 处理 ,而 后 进 
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图 7 局 部 信 噪 比 LLSNR) 和 对 比 度 信 噪 比 (CNR) 感 兴 
趣 区 域 图 

Fig.7 Analysis of LSNR and CNR in the region of 
interest (ROI). The red and blue boxes indicate the 
ROIs for LSNR and CNR analysis, respectively. 


SUCEDE 


行 基于 FBP 的 DBT 重 建 。 重 建 结果 显示 ,在 图 像 细节 
结构 保持 方面 ,本 文 算法 优 于 BM3D 算法 。 

本 文 对 恢复 后 的 投影 数据 采用 了 距离 驱动 反 投影 
(DDBP) 算 法 进行 重建 ,DDBP 算 法 相 比 传统 的 对 整个 
图 像 基 于 路 径 的 反 投影 缩短 了 计算 时 间 忆 ,与 此 同时 ， 
从 精度 的 角度 出 发 ,DDBP 算 法 考虑 了 图 像 像素 和 探测 
器 的 宽度 问题 ,可 以 包 略 将 投影 数据 的 锥 形 束 近似 平行 
束 处 理 时 对 重建 的 影响 。 

总 体 来 说 ,通过 对 不 同 剂量 下 的 乳腺 体 模 数据 的 研 
究 对 比 显示 ,基于 噪声 相关 性 的 惩罚 加 权 最 小 二 乘 算法 
可 以 用 于 低 剂 量 DBT 成 像 投影 数据 降 噪 ,该 应 用 使 低 
剂量 的 乳腺 层 析 图 像 不 再 妥协 于 低 图 像 质量 ,为 乳腺 层 
析 技 术 的 进一步 发 展 提供 了 技术 支持 。 


表 1 不 同 剂量 下 投影 数据 恢复 前 后 的 重建 图 像 对 比 度 信 噪 比 (CNR) 


Tab.1 CNRs calculated from the reconstruction results before and after projection data recovery 


at different doses 


Dose (mA) 
CNR 
40 60 80 100 
Before denoising 0.15 0.21 0.22 0.25 
After denoising 0.63 0.91 1.06 1.29 


表 2 不 同 剂量 下 投影 数据 恢复 前 后 的 重建 图 像 局 部 信 噪 比 LLSNR) 


Tab.2 LSNRs calculated from the reconstruction results before and after projection data 


recovery at different doses 


Dose (mA) 
LSNR 
40 60 80 100 
Before denoising 4.69 537 6.24 6.85 
After denoising 18.55 24.79 30.40 35.49 
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